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ABSTRACT 



We report on the discovery of a molecular cavity in the Norma near arm in 
the general direction of Westerlund 1 (Wdl), but not associated with it. The 
cavity has a mean radial velocity of —91.5 kms~^, which differs by as much as 
~40 kms~^ from the mean radial velocity of the Wdl stars. The cavity is sur- 
rounded by a fragmented molecular shell of an outer diameter of about 100 pc and 
IO^Mq, which is expanding at velocities of 6 to 8 kms~^. The amount of kinetic 
energy involved in the expanding shell is ~ 10^^ erg. Inside this cavity the atomic 
HI gas surface density is also the lowest. Structure of the extended Very High 
Energetic (VHE) 7-ray emission, recently reported by the H.E.S.S. collaboration 



( lOhm et al. 



20091 ) ■ coincides with the cavity. The observed morphology suggests 
that the inner wall of the molecular shell is the zone of the 7-ray emission, and not 
the dense gas surrounding massive stars of Wdl as had been speculated by the 
H.E.S.S. collaboration. A likely candidate responsible for creating the observed 
cavity and the 7-ray emission is the pulsar PSR J1648-4611. 



Subject headings: Gamma Rays: observations — ISM: Molecules — ISM: clouds 
ISM: structure — Open clusters and associations: individual (Westerlund 1) 
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Introduction 



Recently, the H.E.S.S. collaboration carried out a search for Very High Energy 
(VHE, E>100 Gev) 7-ray emission from Galactic young stellar clusters and found positive 
detections in the direction of Westerlund 1 (Wdl) and Westerlund 2 (Wd2) (Ohm et al. 
2009; Ohm09 henceforth). The motivation behind these searches comes from the fact that 
the stellar clusters are potential acceleration sites of VHE particles, since they host a 
variety of energetic sources such as Supernova Remnants (SNRs) and pulsar wind nebula e. 



2003) 



Extended and point like emission was detected surrounding Wd2 (lAharonian et al. 
whereas only extended emission off-centered from the cluster was detected in the c ase of 
Wdl (OhmOQ). The emission surrounding Wd2 was associated with the cluster by 



Dame 



(120071 ). who found structural coincidences between the CO map and the 7-ray source. 
OhmOQ, searched for structural coincidences between the 7-ray and HI maps of Wdl, and 
suggested a possible association of Wdl with the 7-ray emission. One way of confirming 
this association is to compare the observed 7-ray emission structure with the CO gas, which 
is a well-known tracer of the 7- ray emission zones, apart from tracing the spiral arms better 



than HI in the inner disk (e.g. iDame fc Thaddeud (120081 )). 



In order to study the Wdl regio n in CO, we need to establish the velocity of the Wdl 



cluster. iKothes fc Dougherty! (j2007l ). analyzed the neutral gas environment around the 
cluster and concluded that Wdl lies at a radial velocity of —50 km s^^, thus locating it in 
the external part of the Scutum-Crux (SCx) arm. Recent radi al velocity measurements of 



stellar members of Wdl yield a mean velocity of —50 km s ^ (IMengel fc Tacconi-Garman 



2009 



Ritchie et al. 



20091 ). thus confirming the velocity derived from the HI analysis. 



We ha ve searched for molecu lar clouds towards Wdl, in the Columbia- Calan CO survey 



data cube ( Bronfman et al. 



19891 ). This survey covers the entire southern MW and has a 



spatial resolution of 0.125°, a velocity resolution of 1 kms ^ and a sensitivity of 0.1 K. Our 
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analysis covers large angular scale (2° x 2° centered at 1=339.55°, b=— 0.4°), allowing us to 
study the large scale structure of the interstellar medium around the cluster. 

The aim of the present work is to search for the molecular clouds associated with the 
H.E.S.S. 7-ray emission, and to investigate whether Wdl is the source of the observed 7-ray 
emission. With this aim, we first analyzed the morphology and kinematics of the molecular 
gas at the established radial velocity of Wdl (—50 kms~^) in the SCx arm. The resulting 
poor association of the CO structures with the H.E.S.S map, prompted us to carry out the 
analysis of the CO morphology of the Norma arm at —90 kms^^, which is the next major 
arm in the line of sight (LOS) to Wdl. In this latter arm, we discovered a ring, which could 
be interpreted as an expanding shell, with very good structural correspondence with the 
observed 7-ray emission. In this Letter, we compare the CO maps at both velocities with 
the 7-ray observations, and discuss in some detail the CO structure of the Norma arm. 
In §2, the CO molecular emission in the direction of Wdl is presented. The association 
between molecular gas and VHE emission is presented in §3, and finally in §4, we discuss 
the consequences of this correspondence. 



2. CO emission in the general direction of Wdl 



The ^^CO intensity line profiles along the LOS to Wdl in narrow and wide beams are 
shown in Fig. 1. This LOS inte rsects several arms of the MW. We adopted the Galactic 



model of 



Bronfman et al. 



(120001 ) to identify the observed peaks with known spiral arms. The 
3 kpc arm at —120.6 kms^^ with a FWHM of 13.7 kms^^ (distance~7.5 kpc), is identified. 
The peak at —91.5 kms~^ with a FWHM of 8 kms^^ corresponds to the Norma near 
arm (distance=5.5±0.5 kpc), whereas the strongest emission profile in the plot originates 
principally from the SCx arm at —35 kms^^ (distance^ 3.5 kpc). However, velocity 
profiles of Norma far (—50 kms~^, distance^ 12.1 kpc), and Carina arm (—10 kms~^. 




Bronfman et al 



fl2nnnl ) 



Fig. 1. — Top: Reproduction of the longitude- velocity diagram of 
indicating the position of spiral arms along at 1=339.75° (solid line). The flat rotation curve 
are marked with small squares. Bottom: The ^^CO emission profiles in the direction of 
Wdl in a narrow (15' x 15'; solid line), and wide (2° x 2°; dashed line scaled by a factor of 
1/30) beams. The range of radial velocities for 3 of the arms are indicated by the shaded 
boxes. The arrows point to the radial velocities obtained by fitting Gaussian profiles to the 
CO emission peaks relevant to this study. 
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di stance~0.5 kpc), also c ontribute to the observed profile as can be inferred from Fig. 5 



of 



Bronfman et al. 



(|2000[ ). The far side of all arms in the fourth quadrant of the MW, is 
very difficult to detect and/or distinguish, because they are expected to be much fainter 
than the superposed contributions from nearer arms. Hence, we assign the entire intensity 
of a velocity peak to the nearer arm. We identified the peak located at a velocity of 
-48.3 kms-^ (FWHM of 



4 kms~^) with the Wdl cluster, following the HI study by 
Kothes fc Dougherty! ( l2007l ). Given that they used a galactocentric distance to the Sun of 
7.6 kpc, we recalculated the distance to Wdl using the lAU recommended value of 8.5 kpc, 
obtaining a distance of 4.3 kpc (instead of 3.9 kpc). 

The CO emission map centered at the radial velocity of the Norma arm is shown 
in Fig. 2. In the velocity integrated map, an apparent CO ring centered at 1=339.75°, 
b=— 0.4°, is easily identifiable. This ring like structure constitutes a cavity surrounded by a 
molecular gas complex of an outer diameter of 1.5°, and a thickness with an average value 
of about 0.5°. It is part of a fragmented shell seen in projection (see below). 

It can be seen in both position-velocity maps (Fig. 2c and 2d) that there is both 
receding and approaching gas with respect to the systemic velocity of the Norma arm. 
The observed CO emission largely originates from the gas at the arm's systemic velocity. 
The approaching gas is principally concentrated to the north of the ring center, whereas 
the receding gas is located to the south. The approaching gas is not concentrated in a 
single structure, instead the observed emission comes mainly from two regions, one located 
at 1=340.0°, b=0.0° and the other at 1=339.2°, b=0.0°. They have a maximum radial 
velocity of 8 kms^^ with respect to the systemic velocity (—91.5 kms~^) of the Norma 
arm. These components are compact in projection and they merge in velocity smoothly 
with the Norma's systemic gas, as can be seen in fig. 2(d). This leads us to believe, that 
the approaching component belongs to the Norma arm, and not to the 3 kpc far arm. The 
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Fig. 2. — 2D projections of the ^^CO data cube ascribed to the Norma arm (—104 < V^ad < 
—84 kms^^). (a) Velocity integrated intensity map. The lowest gray-scale level corresponds 
to la above the background, with successive levels increasing in steps of 7cr. (b) Contour 
maps showing the receding (red) and approaching (blue) gas with respect to that at rest 
(—94 < Vrad < —90 kms~^; green). The position of Wdl is indicated by an asterisk, (c) 
Latitude- velocity map, and (d) Longitude- velocity map. An area of 1.5° x 1.5° is integrated 
in the projected spatial axis in the last two panels, where we also identify the slices of 
receding, approaching, and rest gas with red, blue and green boxes, respectively. The cloud 
at 1=338.5 and b=0.1 (top-right box in panel c), has been excluded from further analysis. 
The systemic velocity of the Norma arm is marked with an arrow. 
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Fig. 3. — A proposed schematic model based on the CO observation that reproduces the 
integrated maps. Colors identifies the receding, approaching, and rest gas with red, blue 
and green respectively. Z is the depth along the line of sight. The 7-rays emission zone is 
marked with yellow and is discussed in §3). 
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receding component is distributed uniformly along the latitude and longitude range of the 
cloud with a maximum velocities of 6 kms~^ with respect to the arm's systemic velocity. 
Fig. 2b shows these velocity components superposed on the spatial map. All these observed 
characteristics suggest that the molecular ring is expanding at 6 to 8 kms~^. 

In Fig. 3, we present 2D sections of a 3D expanding shell schematic model that best 
explains the observed structures in space and velocity. In Fig. 3 (a),(b),(c) and (c') we 
show the projections in latitude-depth, longitude-depth and latitude-longitude planes, 
respectively. In each panel the approaching and receding gas is shown in blue and red colors 
respectively, whereas the gas at rest is shown in green. The observed spatial distribution of 
gas, as well as the observed kinematics, are well reproduced by our simple schematic model, 
in which the gas principally resides on the surface of an expanding shell, with very little 
gas in the center of the shell. The gas in the shell is not uniformly distributed, instead 
the majority of the gas mass is located in two diagonally opposite directions. From the 
model, we infer that the gas at rest (in green) is also participating in the expansion, but 
it is seen at rest due to projection effects. In Fig. 3 (d) and (e), we show the projections 
in the latitude-velocity and longitude-velocity planes respectively. These figures reproduce 
well the observed structures in Fig. 2 (c),(c') and (d), respectively. 

We estimate the molecular cloud mass content in the shell to b e 10^M(t^. For calcu lating 



1997|), and 



this, we used the ratio N(H2)/Wi2co= 1.56x 10^° cm~^ kms~^ (IHunter et al.l 
with the standard approximation Mcioud={'^) (1.36) A mp N(H2), where A is the projected 
area, m^ the proton mass, N(H2) the molecular column density, and Wco the ^^CO 
integrated intensity, the factor 1.36 accounts for Helium assuming a 10% abundance by 
number (See Murphy & May 1991). The kinetic energy involved in the observed gas mass 
and expansion velocities is about ~ 10^"*^ erg using Efcj„=(l/2) Mc^oMd Vg^.p- The observed 
size of the shell implies that the cavity was produced around 6 Myr ago, if the expansion 
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velocity was uniform over time, or shorter if the shell had a larger expansion velocity in the 
past. 

A map showing the distribution and the kine matics of the atomic HI gas associated 



LunaetaL 



(120091 ). We note that the 



with the molecular cavity is presented in Fig. 4 of 
distribution and kinematics of the atomic HI gas is very similar to that of the CO emission. 
The cavity observed in the molecular gas is seen in the atomic g clS clS cLll HI hole of low 
column density. The lowest HI integrated emission outlines the northern inner boundary of 
the molecular cavity. 



3. The molecular gas associated to the 7-ray emission 

In this section, we compare the morphology of the molecular cavity with the 7-ray 
emission map from Ohm09. The top panel of Fig. 4 shows the CO emission associated with 
the molecular shell in contours superposed on the 7-ray emission map, which has a beam of 
30'. The latter map is dominated by several bright point sources that have been identified 
by Ohm09 as associated with already known sources. The X-ray source GX340-I-0 at 
1=339.56, b=-0.085, the pulsar PSR J1648-4611, and the magnetar CXO J164710.2-455216 
(white asterisk in Fig. 4), are of relevance to the present paper. In addition to these point 
sources, there is extended emission in the shape of a horseshoe, that coincides with the 
projected inner wall of the molecular shell. Especially noticeable is the similarity of the 
structures of the 7-ray emission and the dense molecular gas northward of the magnetar, 
where the projection of the molecular shell is seen as a protrusion inside the cavity. The 
7-ray emission seems to fill the cavity in the figure. Such a configuration can be produced 
in our 3D expanding shell model, if the 7-ray emission comes from the inner surface of the 
molecular shell as can be seen in yellow shades in Fig. 3. Hence, the proposed shell model 
with the 7-ray emission restricted to the inner wall of the shell, can explain not only the 
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Fig. 4. — The extended VHE 7-ray emission detected by the H.E.S.S. collaboration (Ohm09, 
their "preliminary" Fig. 3), where we have overplotted green contours of integrated CO gas 
at rest with respect to Norma arm (yrad=~91 kms^"*^, AV=5 kms""*^ every 3o" from 3a){top), 
and CO gas at rest with respect to Wdl (K,ad=~48 kms~^, AV=20 kms^^ every 10a from 
lOcr) (bottom). Thickness of the line contours are decreased from bright to faint levels to 
remark structures. 
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observed molecular position-velocity map, but also the 7-ray emission structure observed 
by H.E.S.S. 

While the molecular shell of the Norma arm could be traced in many successive channel 
maps. Yet, there is no coordinated structure in anyone the velocity channel maps around 
the systemic velocity of Wdl (~ —50 kms~^). Hence, in order to generate a CO map of 
the Wdl region, we integrated all the channels in a 20 kms~^ velocity range centered at 
—48 kms~^. In the bottom panel of Fig. 4, we show this CO map superposed on the 7-ray 
emission map. There is no coordinated CO structure associated with Wdl, that could be 
compared with the 7-ray emission map. The only prominent structures seen in CO are at 
1=340.25°, b=-0.30° and 1=339°, b=0.5°. Hence, Wdl is not responsible for the observed 
extended 7-ray emission. 



4. Summary and Discussions 

We analyzed the morphology and kinematics of the molecular environment in the 
general direction of Wdl, with the aim of looking for the molecular gas that may be 
associated with the recently reported H.E.S.S. 7-ray emission. We discovered an expanding 
molecular ring at the radial velocity of —91.5 kms~^, with the inner contours of the ring 
coincident with the observed 7-ray emission structure. The observed morphologies of the 
7-ray emission and the CO emission along with their kinematics, can be understood in 
terms of an expanding fragmented shell of molecular gas. Thus, the 2 components are 
physically associated and both lie in the Norma arm. Consequently, the molecular structure 
surrounding the Wdl (—50 kms~^) cluster does not correspond to the observed 7-ray 
emission map. 

The 7-ray emission is produced by accelerated energetic particles when they interact 
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with dense gas. Recently, iFujita et al.l (|2009[ ) proposed a model of generating VHE 7-rays 
by a SNR located in a cavity surrounded by high density molecular gas in order to explain 
the observed photon spectra for the hidden SNR in the open cluster Wd2, and the old-age 
mixed-morphology SNR W28. In this model, the particles are accelerated at the end of 
the Sedov phase, when the SN shock reaches and collides with the surrounding spherical 
high-density molecular gas (shell). This interaction produces 7-ray emission in the inner 
wall of the molecular shell, exactly the configuration that produces our observed CO 
morphology. Hence, the presence of a stellar remnant inside the molecular cavity can give 
rise to the observed 7-ray emission. There are three stellar remnant can didates in the fielc 



of vie w. They are a binary system GX340+0 containing a neutron star (jSchultz fc Wijers 



I993I), a magnetar (CXO J164710.2-455216), and the pulsar PSR J1648-4611, with the 



third one the most likely precursor. 

The first two sources are seen geometrically inside the molecular cavity, but lie at 
distances different from that of the cavity according to their pr esent distanc e estim ations. 
The magnetar is associated with Wdl following the analysis by 



MunoetaL 



mm . based 



on statistical grounds. However, it could as well be located further out, in the molecular 
cavity. For the neutron star in the binary system GX340-I-0, a maximum distance of 
11 ±3 kpc has be en estimated based o n its luminosity, assuming that it is radiating at the 



Eddington limit ( jPenninx et al. 



19931). Yet, it is conceivable that the neutron star is not 
emitting with the limiting luminosity, in which case it could be located at the distance of 
our cavity. 

A third stellar remnant (PSR J1648-4611) is located half a degree South-East of the 
ring center and is seen just at the outer border of the molecular ring (see Fig. 4). A 
distance by dispersion measure yiel ds a value of 5-5.7 kp c, depending on the electron 



density distribution model adopted ( iTorres fc Nuza 



2OO3I ). That is at the same distance 
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of the cavity. This makes this pulsar a strong candidate for producing the 7-ray emission 
at the inner boundaries of the molecular shell. The location of the pulsar off-center of the 
cavity, could be due to a large proper motion of the pulsar, following the SN explosion. It 
requires only a velocity as little as 8 kms~^ for the pulsar to move from the center of the 
shell to its present position, if the explosion occurred 6 Myr ago, and the pulsar is moving 
in the plane of the sky. 

Recently, the LAT instrument on board of th e Fermi satellite detected unpulsed 



em ission from a region coincident with this pulsar (lAdbo et al. 



by 



Fujita et al. 



20091 ). The model proposed 



(I2OO9I ) expects emission in the Fermi band coincident with the 7-ray 
emission, but such emission is not yet detected by Fermi observations. However, the 
observed configuration of a cavity surrounded by a molecular shell, with associated VHE 
7-rays emission, and a bright Fermi source in the vicinity, are common to both, Wd2 and 
the expanding shell discussed in our work. 
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